Orbital-free density-functional theory for metal slabs.

Orbital-Free Density-Functional Theory (OF-DFT) is known to represent a promising alternative to the standard Kohn-Sham (KS) DFT, as it relies on the electron density alone, without the need to calculate all KS single-particle orbitals and energies. Here, we investigate the behavior of the main ingredients of this theory, which are the non-interacting kinetic-energy density (KED) and the Pauli potential, for metal slabs. We derive explicit density functionals for these quantities in the quantum limit where all electrons are in the same slab discrete level of energy, and we present numerical calculations beyond this quantum limit for slabs of various widths. We have found the first explicit KED functional for a realistic many-particle fermionic system, which we prove to be generally valid with no assumption about the KS potential. We also discuss the total non-interacting kinetic energy and the corresponding enhancement factor, which represent basic quantities for the practical implementation of OF-DFT.

Plasmon Modes of Graphene Nanoribbons with Periodic Planar Arrangements.

Among their amazing properties, graphene and related low-dimensional materials show quantized charge-density fluctuations-known as plasmons-when exposed to photons or electrons of suitable energies. Graphene nanoribbons offer an enhanced tunability of these resonant modes, due to their geometrically controllable band gaps. The formidable effort made over recent years in developing graphene-based technologies is however weakened by a lack of predictive modeling approaches that draw upon available ab initio methods. An example of such a framework is presented here, focusing on narrow-width graphene nanoribbons, organized in periodic planar arrays. Time-dependent density-functional calculations reveal unprecedented plasmon modes of different nature at visible to infrared energies. Specifically, semimetallic (zigzag) nanoribbons display an intraband plasmon following the energy-momentum dispersion of a two-dimensional electron gas. Semiconducting (armchair) nanoribbons are instead characterized by two distinct intraband and interband plasmons, whose fascinating interplay is extremely responsive to either injection of charge carriers or increase in electronic temperature. These oscillations share some common trends with recent nanoinfrared imaging of confined edge and surface plasmon modes detected in graphene nanoribbons of 100-500 nm width.

Spherical-shell model for the van der Waals coefficients between fullerenes and/or nearly spherical nanoclusters.

Fullerene molecules such as C(60) are large nearly spherical shells of carbon atoms. Pairs of such molecules have a strong long-range van der Waals attraction that can produce scattering or binding into molecular crystals. A simplified classical-electrodynamics model for a fullerene is a spherical metal shell, with uniform electron density confined between outer and inner radii (just as a simplified model for a nearly spherical metallic nanocluster is a solid metal sphere or filled shell). For the spherical-shell model, the exact dynamic multipole polarizabilities are all known analytically. From them, we can derive exact analytic expressions for the van der Waals coefficients of all orders between two spherical metal shells. The shells can be identical or different, and hollow or filled. To connect the model to a real fullerene, we input the static dipole polarizability, valence electron number and estimated shell thickness t of the real molecule. Our prediction for the leading van der Waals coefficient C(6) between two C(60) molecules ((1.30 ± 0.22) × 10(5) hartree bohr(6)) agrees well with a prediction for the real molecule from time-dependent density functional theory. Our prediction is remarkably insensitive to t. Future work might include the prediction of higher-order (e.g. C(8) and C(10)) coefficients for C(60), applications to other fullerenes or nearly spherical metal clusters, etc. We also make general observations about the van der Waals coefficients.

Van der waals coefficients for nanostructures: fullerenes defy conventional wisdom.

The van der Waals coefficients between quasispherical nanostructures can be modeled accurately and analytically by those of classical solid spheres (for nanoclusters) or spherical shells (for fullerenes) of uniform valence electron density, with the true static dipole polarizability. Here, we derive analytically and confirm numerically from this model the size dependencies of the van der Waals coefficients of all orders, showing, for example, that the asymptotic dependence for C(6) is the expected n(2) for pairs of nanoclusters A(n)-A(n), each containing n atoms, but n(2.75) for pairs of single-walled fullerenes C(n)-C(n). Large fullerenes are argued to have much larger polarizabilities and dispersion coefficients than those predicted by either the standard atom pair-potential model or widely used nonlocal van der Waals correlation energy functionals.

Benchmark quantum Monte Carlo calculations of the ground-state kinetic, interaction and total energy of the three-dimensional electron gas.

We report variational and diffusion quantum Monte Carlo ground-state energies of the three-dimensional electron gas using a model periodic Coulomb interaction and backflow corrections for N = 54, 102, 178, and 226 electrons. We remove finite-size effects by extrapolation and we find lower energies than previously reported. Using the Hellman-Feynman operator sampling method introduced in Gaudoin and Pitarke (2007 Phys. Rev. Lett. 99 126406), we compute accurately, within the fixed-node approximation, the separate kinetic and interaction contributions to the total ground-state energy. The difference between the interaction energies obtained from the original Slater-determinant nodes and the backflow-displaced nodes is found to be considerably larger than the difference between the corresponding kinetic energies.

The Many-Body Exchange-Correlation Hole at Metal Surfaces.

We present a detailed study of the coupling-constant-averaged exchange-correlation hole density at a jellium surface, which we obtain in the random-phase approximation of many-body theory. We report contour plots of the exchange-only and exchange-correlation hole densities, the integration of the exchange-correlation hole density over the surface plane, the on-top correlation hole, and the energy density. We find that the on-top correlation hole is accurately described by local and semilocal density-functional approximations. We also find that for electrons that are localized far outside the surface the main part of the corresponding exchange-correlation hole is localized at the image plane.

Collapse of the electron gas to two dimensions in density functional theory.

Local and semilocal density functional approximations for the exchange-correlation energy fail badly in the zero-thickness limit of a quasi-two-dimensional electron gas, where the density variation is rapid almost everywhere. Here we show that a fully nonlocal fifth-rung functional, the inhomogeneous Singwi-Tosi-Land-Sjölander (STLS) approach, which employs both occupied and unoccupied Kohn-Sham orbitals, recovers the true two-dimensional STLS limit and appears to be remarkably accurate for any thickness of the slab (and thus for the dimensional crossover). We also show that this good behavior is only partly due to the use of the full exact exchange energy.

High-level correlated approach to the jellium surface energy, without uniform-gas input.

We resolve the long-standing controversy over the metal surface energy: Density-functional methods that require uniform-electron-gas input agree with each other, but not with high-level correlated calculations such as Fermi hypernetted chain and diffusion Monte Carlo calculations that predict the uniform-gas correlation energy. Here we apply the inhomogeneous Singwi-Tosi-Land-Sjölander method, and find that the density functionals are indeed reliable (because the surface energy is bulklike). Our work also vindicates the use of uniform-gas-based nonlocal kernels in time-dependent density-functional theory.

Hellman-Feynman operator sampling in diffusion Monte Carlo calculations.

Diffusion Monte Carlo (DMC) calculations typically yield highly accurate results in solid-state and quantum-chemical calculations. However, operators that do not commute with the Hamiltonian are at best sampled correctly up to second order in the error of the underlying trial wave function once simple corrections have been applied. This error is of the same order as that for the energy in variational calculations. Operators that suffer from these problems include potential energies and the density. This Letter presents a new method, based on the Hellman-Feynman theorem, for the correct DMC sampling of all operators diagonal in real space. Our method is easy to implement in any standard DMC code.

Role of surface plasmons in the decay of image-potential States on silver surfaces.

The combined effect of single-particle and collective surface excitations in the decay of image-potential states on Ag surfaces is investigated, and the origin of the long-standing discrepancy between experimental measurements and previous theoretical predictions for the lifetime of these states is elucidated. Although surface-plasmon excitation had been expected to reduce the image-state lifetime, we demonstrate that the subtle combination of the spatial variation of s-d polarization in Ag and the characteristic nonlocality of many-electron interactions near the surface yields surprisingly long image-state lifetimes, in agreement with experiment.

Hole dynamics in noble metals.

Hole dynamics in noble metals (Cu and Au) is investigated by means of first-principles many-body calculations. While holes in a free-electron gas are known to live shorter than electrons with the same excitation energy, our results indicate that d holes in noble metals exhibit longer inelastic lifetimes than excited sp electrons, in agreement with experiment. The density of states available for d-hole decay is larger than that for the decay of excited electrons; however, the small overlap between d and sp states below the Fermi level increases the d-hole lifetime. The impact of d-hole dynamics on electron-hole correlation effects is also addressed.